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SYNOPSIS 


la this wr] angular distributions a ul excitation curves 

IQ Q 

of (d,p) reactio i on C and Be have been obtained at the 

-.ob-coulomb v 'onDarc'in g energy in the ra^ge of 500 to 1800 KeV 

using the 2 MeV ? an-ae-Graaf f accelerator at I I T Fanpur 

12 

Strong resonances were observed m the case of C and compa- 
ratively weal er resonances m the case of ^Be m the excita 
tion curves The angular distributions in both cases however 
showed a peak at about 35° M typical of direct stripping mode 
At resonant energies the main effect was in the enhancement 
of the backward peak Some data on both these reactions was 

previously available However, the analysis was done either 

m terms of compound nucleus process alone or the stripping 
mode alone In the case of compound nucleus mechanism, it was 
necessary to assume contribution from rather high partial waves 
and m the case of stripping mechanism, it was necessary to 
assume optical model parameters of DWBA calculations which 
were not consistent with those used at highei energies where 
the angular distributions show better stripping patterns It 
was therefore, decided to analyse the data m terms of coherent 
contribution of the reasonant and non-resonant amplitudes and 
then compare the non-resonant cross-section with the DWBA cal- 
culations with optical model parameters consistent with those 
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used at higher energies For a proper consideration of the 
interference effects it was necessary to ia e the data at 
reruired energies and angles although so c oils was availa- 
ble from previous work 

The amplitude was assumed to be of the i->rm 
^ _ ^Resonant ^Non-resonant 

m the vicinity of resonance peal there ^ es01 ant included con- 
tribution from tne main peak and one or two neighbouring peaks 
It was assumed that the energy dependence of the resonant 
amplitudes was given by Breit-Uigner tern with phases of the 
neighboars relative to the mam resonance The airect ampli- 
tude was considered relatively ansensitive to energy over the 
region of the resonance and was assumed to be at the most 
linearly dependent on it with the undetermined phase relative 
to the main resonance This assumption about the total ampli- 
tude leaf s to an exnressioh for the total cross-section with a 
term corresponding to the non-resonant cross-section, and 
others depending upon energy m the specified wa^ The expe- 
rimental cross-section m the neighbourhood of resonances m 
the excitation curve was then subjected co a least squares fit 
using the expression for cross-section with strengths of various 
terms resonant energies and widths, as parameters The least 
squares fit provides a value of the non-resonant cross-section 
at each energy of the excitation curve obtained at each angle 



It is then possible to reconstruct tne angular distribution 
due to the non-resonant process alone 

T hese extracted non-resonant angular distributions were 
compared wifi a D B \ calculations of the stripping process 
made vith a program DWUCK available at the rhabha Atomic 
Research Centre It was found that there were satisfactory 
agreemert with the shape and absolute cross-section between 
the two 

II is therefore concluded that while the contributions 

of various processes m determining the cross-section are to 

be determined by coupled ch-'ir^l calculations the simpler 

model of adding resonant and non-resona^t amplitudes suggested 

by physical considerations m the classical limit seems to 

yield satisfactory results l ith this approach the contnbu- 

12 

tion of the direct reaction m the case of C was found to be 

o 

about 20 / and that m the case of ""Be about 90 / In view of 
tr is the assignment of spins and parities made previously by 
neglecting the non-resonant part needs to be re-examined An 
attempt to fit the extracted direct cross-section with Coulomb 
distortions alone did not provide satisfactory fits mspite of 
the low bombarding energy Nuclear distortion is therefore, 
found important even at these energies The finite range and 
non-local corrections the DWBA calculations did not produce 
any major change m the shape however the absolute cross- 
section reduced by about 10 / with these corrections applied 
to the zero range approximation 


CHAPTER I 


INTRODUCTION 

The striping feature that was noted m the early stages 
of nuclear reaction worp was the appearance of narrow resonances 
m the cross-section as a function of bombarding energy This 
led to the suggestion of the compound -nucleus mode of nuclear 
reactions by N Bohr 1 with compound -nucleus life time of the 
order of 10 10 sec which was much larger than the typi- 

cal transit time of the bombarding particle across the nuclear 
diameter The theory of nuclear reactions on the basis of com- 
pound-nucleus formation mechanism was subsequently developed 

2 

by many authors and was successfully applied to determine the 
properties of many compound nucleus levels 

In the case of nuclear reactions proceeding via compound 
nucleus formation mechanism it is expected that excitation fun- 
ction will show resonances whose shapes are given by the Breit 
Wigner formula Since the formation and decay of the compound- 
nucleus are independent of each other the same reaction products 
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would be obtained even if a particular level of the compound- 
nucleus is formed by two different initial channels The angu- 
lar distribution of the reaction products proceeding through a 
single level is expected to be symmetric about 90° m the centre 
of mass system Since the shape of the angular distribution is 
characteristic of the partial waves m the outgoing channel it 
is possible to make an inference about the spin of the compound- 
nucleus level If at relevant excitation energy two or more 
levels of the compound nucleus of different panties contribute 
the angular distributions are not expected to be symmetric 
and the shape of the angular distribution is expected to vary 
noticeably with the bombarding energy The width and the den- 
sity of the compound -nucleus level are expected to increase with 
excitation energy and when a large number of levels contribute 
it is expected that, due to essentially random phase of contribu- 
tion the angular distribution would again be symmetric and 
weakly dependent on angle 

Even before the suggestion of compound -nucleus process 

1 3 

by Bohr Oppenheimer and Phillips had proposed a radically 
different model to account for the higher probability of (d p) 
reactions as compared to (d n) reactions They pointed out 
that at the low bombarding energy which was used at that time 
the proton in the deuteron would not cross the Coulomb-barrier of 
tne target while the neutron would get close enough to get cap- 
tured It was thus possible to explain the preferential proton 


3 


emission which is lust the opposite of the predictions of com- 
pound-nucleus theory The full implication of this model was 
not clear at that time since at very low deuteron energies the 
angular distributions were isotropic as only s-wave interactions 
took place In the fifties (d p) and (d, n) reactions were 
studied at higher energies and they were found to be equally 
probable showing that Coulomb force p] ays the secondary role 
However the most interesting feature was the angular distribu- 
tion The cross-section peaked either in the forward direction 
or at some small angle and decreased to small values at back- 
ward angles These features were independent of bombarding ener- 
gies, and peak position depended upon the target used These 

4 

features were explained first by Serber and later more accura- 

5 

tely by Butler The explanations were based on the single step 
model of the reaction in which the deuteron is stripped of one 
of its nucleons, due to the low binding energy of the deuteron 
during the close passage of the captured nucleon relative to the 
tafget nucleus The other nucleon can still have a relatively 
large impact parameter m view of the large deuteron radius 
In the stripping process therefore the deuteron as a whole and 
the target do not combine to form a compound nucleus Butler 
showed that the angular distribution is essentially of the form 
of the square of the spherical Bessel function 5 6 J^(QR) of order 
1 where 1 is the orbital angular momentum transfered to the 
target nucleus by the captured nucleon Q is the linear momentum 
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transferred and R is the nuclear radius The study of deuteron 
stripping reactions therefore provided a new technique of nuclear 
spectroscopy to study the properties of the low lying levels of 
the residual nucleus, while the reactions proceeding through 
compound nucleus mechanisms led to the study of highly excited 
levels of the compound nucleus The inverse (p, d) and (n, d) 
reactions are called the pick-up reactions m which the inci- 
dent nucleon picks up an appropriate peripheral nucleon to form 
the emergent deuteron without forming the compound nucleus The 
Bulter theory is directly applicable to such reactions and the 
peak position in the angular distribution now depends upon the 
angular momentum of the picked nucleon in the target nucleus 
The pick-up reactions thus provide information about the ground 
state of the target nucleus 

In a rearrangement collision the amplitude of the reaction 
7 

a-# Pis given by 


Ace - / V e V B ¥dT <1 1) 

where is the internal wave function of the reaction products 
m channel g the final state interaction and ¥the total 

wave function In the (d p) reaction Vg = V + where 

Vp n is the proton -neutron interaction and is the sum of 

the interactions of the proton with the nucleons of the target 
In the pure stripping model it is assumed that the interaction 
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, ga 

V can be negletced The Butler treatment and the other 

8b 

alternative treatments that appeared soon afterwards are essen- 
tially equivalent to making a plane wave Bom approximation 

(PWBA) m which the total wave function 'Pis replaced by plane 
lF r 

wave e ■ o m the incident channel The method used by 
8 b. 

Butler involves wave function matching at the nuclear boundary 
and it leads” to same result as the PWBA with the difference that 
the radial integral for amplitude is cut-off at the nuclear radi- 
us R and does not extend to the origin 

In spite of the success of the PWBA theory for deuteron 

stripping reactions, the theory does not yield a particularly 

detailed account of the angular distributions For medium and 

heavy nuclei the agreement between the PWBA theory and the expe- 

g 

rimental angular distributions is poor Even for light nuclei 
there are cases m which the assignment of the orbital angular 
momentum is ambiguous ^ Further PWBA gives a value of the 
absolute cross-section which is several times the measured cross- 
section 

The exact expression (1 1) for the amplitude of the 
reaction can also be put m the form 

= /4~ )( ^>fp ( v e -Vp) * a T (l 2 ) 

where Vg is a model potential and Xg(r^ ) represents the relative 
motion m the model potential The PWBA is improved upon 
by using the Distorted Wave Bom approximation (DWBA) theory 
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The total wave function ¥ m expression (12) is replaced 
not by the plane unscattered wave but by the scattered wave 
X„ (t" ) m the incident channel The scattering state wave funC' 
tion is generated from the optical model where the optical 
model parameters are determined from the available scattering 
data m the incident channel The model potential is chosen 
equal to the optical potential which gives correct scattering 
m channel g Thus ) is the scattering-state wave function 

m channel 6 This wave function is generated from the optical 
model potentials which determine the optical model parameters 
for a given system As mentioned before for the (d, p) reac- 
tion V 6 =V pn + V p(J therefore 

V S--^6 = v pn + ( V" U 3) 


It is assumed that is a close approximation of V 


pN ' 


which is the total interaction the proton experiences due to 
the nucleons m the target nucleus Thus m the DWBA theory the 
plane wave is replaced by the distorted wave which is a better 
approximation to the total wave function and instead of neglec- 
ting the interaction of the outgoing nucleon with the target it 
is assumed that V ^ is equal to the optical interaction to a good 
approximation 


The optical model potential, which is m general required 


m DWBA treatment includes both Coulomb and nuclear interactions 
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The inclusion of nuclear interaction m both the incoming and 
outgoing channel represents a departure from the pure stripping 
mechanism as envisaged originally by Butler It implies that "the 
deuteron as a whole does experience nuclear potential Similarly 
the outgoing proton does not entirely pass through the region 
exterior to the nucleus However the optical potential is 
adjusted to give the correct energy -averaged scattering 

amplitude and does not reproduce the fluctuations m the elas- 
tic channels The potential Vg is not identical with V „ 

DWBA treatment/ therefore/ does not take into account the possi- 
ble compound nucleus formation and subsequent contribution to 
the reaction leading to narrow resonances m the cross-section 
DWBA however provides a more accurate description of the dyna- 
mical behaviour of the incoming deuteron and outgoing proton 
and by doing this provides a sounder basis for using proton- 
neutron interaction as the operator for the matrix element m the 
amplitude (1 2) 

Generally speaking DWBA calculations require extensive 

use of computers Several computer programmes for such calcu- 

13 

lations are now available Many programmes use a zero-range 

approximation for the neutron-proton interaction since it re- 
duces the dimensionality of integral (1 2) Calculations 
which take into account the finite-range of V are lengthier 

ir Li 

Programmes which include such finite range correction as well 
as the spin-orbit term m the optical potential are also 
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available Computations are also made with non-local poten- 
tials There is usually a provision to use a cut-off radius 
m integral (12) as an option The DWBA calculations usually 
give a much more satisfactory account m the entire range of 
angular distribution as compared to PWBA calculations Further 
the absolute cross-sections predicted on the basis of DWBA cal- 
culations are m far better agreement with the experimental 
results 

For bombarding energies which are lower than the Coulomb 
barrier the (d p) cross-sections are expected to be small and 
the Coulomb distortion is expected to be the predominant effect 14 
Secondly the relative contribution to the reaction due to compound 
nucleus formation is expected to be larger This is due to the 
fact that the relative strength of the low partial waves m the 
incident wave increases as energy is decreased and the compound 
nucleus formation is more probable with such low partial waves 
with correspondingly low centrifugal barrier On the other hand, 
when compound nucleus is formed at high bombarding energy it is 
highly excited and therefore has several open channels The 
fractional contribution of the decaying compound nucleus m the 
outgoing proton channel is therefore reduced At low bombarding 
energies, therefore, it is expected that the cross-section would 
be low, distortion effects will be primarily due to Coulomb field 
and the compound-nucleus process will make a significant contri- 
bution to (d,p) reaction 
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The Q value of (d p) reaction is an important conside- 

15a 

ration particularly at low-nombarding energies When the Q 

value of the reaction is large the momentum of the outgoing 
proton is large To have this large momentum it is expected that 
the proton is close to the neutron m the deuteron when the 
neutron is stripped Therefore, such protons pass close to the 
nuclear region and experience marked nuclear distortion It is 
expected, therefore that m the high 0 reaction the plane wave 
approximation will be poor This has been observed for several 
cases 15bc Tne (d,p) reaction on 12 C and 9 Be have been 
studied at bombarding energies xn range 0 8 - 2 0 MeV In these 
reactions Q values are high Therefore both compound -nucleus 
effects and high distortion effects m (d p) reactions are 

expected 

When both compound -nucleus and direct amplitudes are 

present it is expected that they will add coherently While 

. -21 

the duration of a direct event should be of the order of lO 
seconds and that of the compound -nucleus some 10 4 times greater 
the duration of the wave packet which is incident from the acce- 
lerator is of the order of 10 seconds The difference in 
delay between the direct and compound-nucleus mechanism is 
therefore insignificant The well known experimental evidence 
of the existence of interference is that between the potential 
scattering and resonant scattering of neutrons The time 
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intervals involved m direct interactions are the same as those 
nvolved in potential scattering Uhen the interfering ampli- 
tudes are feu m number the interference effects will be noti- 
ceable m the cross-section measurements When the number of 
interfering amplitudes is large the interference terms may tend 
to car cel 

In the presence of an isolated resonance m the (d,p) 

o *t” 

reaction one nay assume the presence of two amplitudes q 
and q adding coherently If the resonance ^s not 

well isolated one may have to add the amplitudes due to neigh- 
bouring resonances In the statistical region where a large 
number of resonances contribute, this approach is not possible 
In the analysis of excitation function m the presence of one 

or more resonances, the resonant amplitudes may be assumed to 

2e 

vary according to the Breit-Wigner formula while the Non- 
Resonant (NR) amplitude may be assumed essentially constant 
over the width of the resonance It is then possible to extract 
the NR cross-section from the experimental excitation curves by 
fitting them with the algebraic function written down on the 
basis of above assumptions The extracted NR cross-section can 
then be analysed on the basis of DWBA 

In the present work we have studied deuteron induced 
reaction at low bombarding energies on light targets with the 
2MeV Van de Graaff accelerator facility available at the Indian 
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Institute of Technology, Kanpur 12 c and 9 Be were used as 

targets Previous data on carbon-12 (d p) reaction m this 

17a 

energy range have been published m 1960 which showed pre- 
sence of resonances and the data was anlysed by the authors 
m the frame-work of compound -nucleus process alone although 
the angular distributions did retain some feature of the main 
stripping peak which shows up more prominently at 4 MeV 16/18 
A preliminary calculation of pure stripping at this bombarding 
energy showed that the stripping cross-section may not be negli- 
ble It was, therefore decided to study the contribution of 
the direct mode relative to the compound-nucleus mode, by consi- 
dering coherent addition of amplitudes For this it was nece- 
ssary to check previous data and to take additional data to 

study the interference effects Some data on (d,p) reaction on 

9 1 q 

Be m this energy range was also available which was analysed 

by the authors m terms of compound -nucleus process alone or 

20 

the direct reaction process alone using optical parameters 
not consistent with those used at higher energies 

12 

In the present vorK it was found that in C(d,p) reac- 
tion the direct mechanism makes about 20/ contribution as compa- 
red to the figure of less than 2/ mentioned m the earlier 
work 1 The extracted MR cross-sections show strong distor- 
tions as expected m relatively high Q value reactions 
The extracted cross-section can not be fitted by pure Coulomb 
distortions. Reasonably good fits are obtained with optical 
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potentials with parameters consistent with those used at 
higher energies In the case of beryllium-9 the compound- 
nucleus contribution is found to be small The extracted 
cross-sections can again be fitted reasonably well with opti- 
cal model parameters consistent with those used at higher 
energies 
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NUCLEAR REACTIONS 


Introduction 

In a nuclear reaction a target nucleus 'A' with n^ 
number of nuclei is bombarded with a projectile 'a' consis- 
ting of n number of nuclei and the reaction products are 

3, 

detected m the asymptotic region The reaction 

a + A 3*" B + b * (2 1) 

is described by denoting the initial and final channels by 
a and 6 respectively Besides the identification of initial 
and final products other relevant quantum numbers are denoted 
collectively by symbols l and j respectively, so that the reac- 
tion is represented by (a,i) — ► (ft,j) when necessary 

The total system thus consists of n,+ n = N particles 

n 3. 

with the Hamiltonian H and a wave function ^ in practice the 

state is a non- stationary state, to be represented by colli- 

21 

sion of wave packets but it is generally represented as a 
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stationary state satisfying 
H X = E X 


(2 2 ) 


The stationary state wave- function ¥ can he expanded 
m terns of a suitable basis such as independent particle 
model (IPM) states of N particles It may be noted that for 
sufficiently high excitation energies the finite well rpM pro- 
vides both bound and unbound states The excitation energy- 
can either be distributed over many excited single particle 
bound states or it may be m the form of one or few particles 
excited into the single particle contimum The bound states 
of IPM therefore provides a discrete spectrum m the same 
region where the unbound states provide continuum The unbound 
states represent the flux of particles m the asymptotic region 
In nuclear reactions the energy of the total system is suffi- 
ciently high to provide unbound components at least m the 
initial channel The wave-function I could then be written in 
the form 


X 


¥ + 
round 


X 


unbound 


(2 3 ) 


^here Y, , is a summation over the 
bound 

j is a summation over unbound 
w unbound 

is often expressed as 


bound basis states and 
basis states l^^a 


X 

unbound 



(2 4 ) 
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where $^is the product <^4^ of the internal wave functions of 

clusters b and P and £ (r„ ) describes their relative notion 

B 3 

While this description of ! , , could be used for large 

unbound 

values of r^ it is sometimes used for snail values of the 
order of nuclear radius as a model wave function It is assumed 
that m the asymptotic region where the boxind states do not 
contribute the total wave function has the form 



lK tf ? cx 

$ e 
a 



U 3 

2 % A a& 

2irn 


lK 3 r 3 

e 


r 3 


(2 5) 


This expression describes an incident plane wave m 
channel ct and outgoing spnerical waves m various channels 3 
Since orly the unbound states contribute m the asymptotic 
region the expans lor (2 4) must have the asymptotic form (2 5) 
The quantity A a g (K^ lg) is called the amplitude of the reaction 
a -* Band the cross-section is proportional to the square of 
this amplitude Measurement of cross-section therefore provi- 
des a sampling of the unbound component of the total wave- 
function 


Ht , and ¥ , , , however are two components 

bound unbound 

of the total wave function and these are directly coupled with 
each other through the Schrodmger Eq (2 2) Thus, if we 
introduce projection operators P and Q with the property 


P ^ ^unbound ^ ^ ^bound 


(2 6 ) 
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5g(rg) is empirically determined to fit the observed data m 
the corresponding scattering channel alone It is m general 
not possible to find model optical potentials which will re- 
produce narrow fluctuations typical of a system with many 
degrees of freedom 


The approach of the compound nucleus theory is to assume 

that while ¥ , , is the dominant term of ¥ m the region 

unbound ^ 

starting just outside the nuclear volume ^ound lS ^ om;L “ 

nant term inside the nuclear volume The coupling between the 
two is then taben into account by the physical requirement of 
matching the two wave functions at the nuclear surface The 
fluctuations can then be understood m terms of rapidly chan- 
ging boundary conditions as determined by ¥k 0unc j An assum P"“ 
tion of smooth behaviour of the logarithmic derivative on the 
nuclear surface then leads to Breit-Wigner formula for the 
resonance The resonance is attributed to a predominant exci- 
tation of a bound state component 


Both the treatments involve approximations to the actual 
problem end introduction of models within these approxima- 
tions The support for introduction of such simplification 
however, is provided by the experimental data which m diffe- 
rent suitable circumstances correspond to good approximations 
to the two different treatments Further, the physical consi- 
derations at least m the classical limit lend support to the 
introduction of these simplifications 
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In a semiclassical approach the nuclear reaction can 
be considered m terns of the development of rearrangement 
collision m time Thus when a deuteron is incident on the 
target there is a possibility of the neutron being captured 
and the proton released Alternately/ both the nucleons may 
enter the target form a compound assembly where they excite 
one, few or many nucleons and eventual] y a proton is released 
again The life time of the compound state may be expected to 
be small if fewer nucleons get excited by the incoming deuteron 
while that of a state which is complicated may be expected to 
be relatively long One would therefore expect resonances with 
varying widths on the basis of uncertainty principle In the 
event that the probability of the simplest process, where the 
proton does not enter the nucleus at all is large compared to 
others, a treatment purely on the basis of direct reactions 
might give satisfactory results In general however both pro- 
cesses will contribute to the observed proton yield This is 
the basis for adding the amplitudes as calculated on the basis 
of two respective mechanisms In general the two amplitudes 
have to be adaed coherently when the wave-packet representing 
the incident particle is much longer than the expected time 
delay caused by the direct and compound nucleus reaction 

mechanisms A resonance of width aE is observed with the 

K. 

beam of energy resolution Mb when 


AEg is smaller than 
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The corresponding uncertainty At_. m the arrival of the beam 

.D 

particle is larger than the time delay At which is of the 
order of nuclear life time Hence the dimension of the beam 
packet is lar^e compared to the distortion produced at its edges 
by the finite life time of the compound nucleus It is expec- 
ted therefore that the compound -nucleus amplitude will add 
coherently and produce interference effects The experimental 
evidence of this situation is provided by elastic scattering 
experiments where one finds interference effects between poten- 
tial and resonant scattering m the excitation curve 

Direct Reactions 

In the direct reaction theory the Hamiltonian of the 
system is considered m the form’*’ 1 


H 

a 

+ 

h a 

+ 

T bA + V aA 

(2 

9a) 

H a 

+ 

T a 

+ 

v a 

(2 

9b) 

H e 

+ 


+ 

Vg 

(2 

9c) 


where H a is the sum of the internal Hamiltonians of the two 
clusters a and A T a the relative kinetic energy of the two 
clusters and is the sum of the potential energies of all 
particles m the cluster a due to each particle m cluster A 
The form (2 9) of the Hamiltonian is not suitable for conside- 
ration of bound wave functions but is suitable for the 
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consideration of unbound wave functions Let$g be the eigen 
function of Hg so that 


H e = e 3 * 6 (2 10 > 

The Schrodmger equation for the total wave function is 

H y (+) = E y (+) (? 11) 


taking the scalar product with we have 
<*g H W ( + ) ) = E (* g , Y (+) ) 


(2 12 ) 


where the integration is only over the internal variables of 
$0 Substituting Eg (2 9c) for H into (2 12) we get 


(Tg+e g -E) ($g * ( + ) ) = -(* p Vg W ( + ) ) (2 13) 

(+) 

Since the scalar procuct ( V ) involves integration over 
only the internal variables of , it defines a function of 

P 

r*g , which is denoted by 

<v s K e <? b > 

Therefore we have 


(t 6 + V E > S “ V 8 ' f(+>) 


(2 14) 


This equation has solution 


Sf}(%> = e 


lK a r a 


2 ,h 2 


iK $ (r B -r B ) () __ 

e~ — — V g ¥ (+, )drJ 


r e " r e 


(2 15) 
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where we have 


E -e 


( h Kg) 

e = 27T~ 


Taking the asymptotic term of eqn (2 15) , the amplitude 


A a6 = - 




2 TTh 


( e 


iKg r 6 


V g ) 


(2 16) 


The formula represented by Eq (2 16) is exact It can 
be recast into a different form which is useful for the dis- 
cussion of applications We add the term Vg ( $g ) to both 

sides of Eq (2 13), where Vg is an arbitrary function of r p 
This leads to the relation 


(T p + V 0 + Eg - E) Sg = -(* g (Vg-Vg)<r (+) ) (2 17) 

The solution of this eqn is obtained by using the Green's 
function technique as above This leads to the following 
expression for tne amplitude 


K„) = - -^-2 ( V B ( ) (V e - V V<+>) 


(2 18) 


2 TT h 


where is a- solution of the homoqeneous Eq 

P 


(Tg + Vg + Eg -E) Xg = 


0 


(2 19) 


Both forms given by Eq (2 17) and (2 19) are equivalent For 
applying either of the two forms a knowledge of the total 
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(+) 

wavefunction V is needed In the direct reaction theory 

the bound state component of is neglected From amongst 

the unbound components often only the component m the incident 
scattering channel is retained Use of more than one unbound 
component leads to a more complicated coupled channel calcula- 
tion The calculation further requires the use of the channel 
potential Vg Simplifying assumptions are usually made for 
calculations 

In tie Plane Wave Bom Approximation (PWBA) the total 

wave function is replaced only by the incident plane wave 

iK(x r a 

<S> a e in the incident channel a The scattered part of 

the incident wave os therefore further omitted In the case 
of (d p) reactions the approximation to Vg is made as follows 
We have 


Vo 


V 


pn 


+ V 


pN 


where V is the proton neutron interaction and V is the 
pn pN 

total interaction of the outgoing proton with all the nucleons 
of the target In the pure stripping mechanism model of Butler 
it is assumed that the proton does not enter the nucleus at all 
It is therefore assumed that the contribution due to the inter- 
action V could be neglected m evaluating the integral (2 16) 
In this approximation the amplitude becomes 
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A a g - - 


2 Fh 


- / 4 1 e 
2 ; g 


-ik fi 


ik„ 


V e 

np a 


dx 



f <4> 


B 


r A 





’d 


dr dr 
n p 


(2 20 ) 


where 4> and ^ are expressed m terms of the proton and the 

deuteron coordinates resoectively The integration in the 

scalar product <<(»_ | 4»_ > is over the internal coordinates of the 

o A 

target nucleus A and it has the form 



s ru (r _ T ) 
1 nN 


(2 20a) 


where s is essentially the coefficient of fractional parentage 
and (r ) is a normalised function of the neutron coordinate 
typical of the angular momentum 1 carried by the neutron bound 
to the target N 

Further the interaction V m Ea (2 20) can be written 

m terms of deuteron binding energy by writing the Schrodmger 

equation for the deuteron Finally changing the variables from 

r / r to the relative coordinates r and r XT we get 
n p pn nN 3 


O V „ -iq r 

2 v pN r f a /- \ _ npy 


A= s ( t(J + q dr np ^(r^) e 


x {lSr nN n * (3 W e 


iQ r 


nN 


(2 21 ) 
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where q and Q are essentially the change of momenta of the 
proton and the nucleus with some scale factors In order 
to evaluate the integrals m (2 21) the exponential factors 
are expanded m terms of spherical Bessel functions as follows 

e 1 ^ r = 4ir l i X (qr)Y* (OjY* (<J ) (2 22) 

X \i 

Since involves spherical harmonics of order zero and ni 
that of order 1, the integration over angles leads to the 
appearance of Bessel functions and respectively For 
the evaluation of integral over often a cut off radius R of 

the order of nuclear radius is used The angular dependence 
of the amplitude is then primarily given by 

A * J ]_ (QR) (2 23) 

The angular distribution of the proton is therefore 
determined by the value of 1 and the momentum transfer Q 
which depends upon the angle of emission The radius R is 
usually used as a parameter within limits The Bessel func- 
tion has an oscillatory behaviour with rapid damping The main 
peak of | J 1 (QR) ] is compared with the experimental main peak 
m (d p) reactions The peak position m many (d p) reactions 
is predicted to a good approximation by this theory but not 
m all cases The agreement away from the peak is however 
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not found satisfactory For prediction of absolute cross- 
sections estimate of the coefficient s is obtained from 
structure considerata on ana using this the predicted cross- 
section is found to be higher than the observed cross-section 
by a large factor 

In the Distorted Wave Bom Approximation (DWBA) the 
total wave function + ^ is replaced by the incident channel 
wave function including the scattered part To obtain this 
incident channel wave- function , scattering data m the 
incident channel is used The parameters of the optical model 
potential are first determined using scattering data and 
these are then used to calculate the incident channel wave- 
function In the expression (2 18) for the amplitude the func- 
tion Vg is any arbitrary function of the channel radius r & and 
the function is then an eigen function of the two body prob- 
lem with the potential Vg In the DWBA the function Vg is 
chosen equal to optical model potential m channel g which 
leads to the observed scattering data m the channel g For 
the (d,p) reactions we then have 

= V + v pN - V opt (2 24) 


It is then assumed that the last two terms cancel each other 
to a good approximation, so that the remaining interaction is 


then V. 


pn 
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The interaction V . T is the sum of all interactions 

pN 

that the individual nucleons of the target exert on the out- 
going proton It is not clear that this sum can be represen- 
ted by a function of the channel radius alone Such reduction 
of a many body interaction to an equivalent two body inter- 
action is an approximation However, the reduction of the 
operator (Vg-Vg) m (2 18) to V pn by a proper choice of V p is 
a much better approximation than neglecting V as is done m 
the FWBA Thus the DT T BA method involves a better approxima- 
tion to i/ + }as well as to the operator m the matrix element 
The amplitude m Eg (2 2o) is then givren by 


k aB 




2irlT 


1 % *3 


(-)* (+) 

K J V $ X dT 

pn a *a 


(2 25) 


where the functions X^are the solutions of 


- 2 , 2 
• ft k 

" 2 U« 


T - V 6 } ^ 0 


(2 26) 


with the asymptotic forms 




lk r ^ , 

■* e + f (e. 


ikr 


— — * e 1 ^ r + f ( it — 0 , tt + 4>) 


-ikr 


(2 27) 


The optica^ pot ntials ^ g arc of the form 


v e = u c + u h 
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where U is the Coulomb potential m the channel and U, T is 
the optical potential This potential is of the form 

and various forms for the real and the imaginary part are 
commonly used 

The amplitude (2 25) can be written m the form 

A o6 = -J ~ / [dr x*" 5 <s|V | a » (2 28) 

2 irn 

where J is the Tacobian of transformation to the relative co- 
ordinates and <£ | V ( a> includes integration over the inter- 
nal co-ordirates of the target 


Iv | a> 

pn 




V 

pn 




a 


n(r T ) v 4 (r ) 

nN pn a np 


This factor which depends upon r_^ and r can be expressed 
m terms of the channel radii r 0 and rg The integral (2 28) 
is then a six dimensional integral Often the interaction V 

pn 

is assumed to be of zero range and is replaced by a delta 

function D 6 (r ) This effectively leads to the 
o np 
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simplification/ 

- - M B - 

r n = r a = ** 

So that the integral (2 28) reduces to three dimensions 

In evaluating the integral (2 28) the matrix element 
<$|Vp n |a> is expanded into a series of multipoles each of 
wmch corresponds to transfer to rhe target nucleus a definite 
angular momentum j composed of an orbital part 1 and a spin 
part S We have then 


J = 1 + s = J B - J A 

and in the absence of the spin-flip we have 
? = ? a - \ 

Often only one value of j , 1 and S is important or allowed m 
a transition The distorted waves x are also expanded m terms 
of partial waves, 


X = 


4 7T 


L, Y M L (r) Y m L (k) ^ X 


kr LM M 


M 


we then have 


1 “ L a “ L e 


The cross-section which is proportional to the square of the 
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amplitude can then be written m the form 
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do _ 2 j b +1 j 


Ks/ 


a ls J 
DW 


d^ 2J a +1 i sj 1 Qx10 4 (2 j+1 ) 


(2 29a) 


where 


a lsj 
DW 


v 

4ir r a 


E a E b y b 


M 


: B 


M 


'A 


1 0x10 _ n \ r 
2S +1 ( 21+1 ^ 

3 . 


™a m b 

3 lsj 


(2 29b) 

where B, is the measure of strength and is proportional to the 
is J 

, mnunv 

spectroscopic factor s, The quantity S a D is given by 

s J l s J 


, mrn a m b 
’lsj 




L -I , -1 
a b 


T L L 
a b a o 


(L a s a° m al J a m a ) <L b s b° m bl J bV 


( J b J rn a -rnm i J a m a^ ( 2L b +1 ) ( Ij a 1 00 I L_o) v&s +1) (2j+l) 


\/(2 J-j^+1 ) (2L a +l ) 


P-b s b J b "7 

1 s j 
L a s a ^a 


d <V 

m b m a" m 


("ir 1 dr a ^ taj'V x L a (k a r a >) 


B 


(2 30) 


J b ^ 


where m = rru - nx + m, - m and d ~ (k, ) is the rotational 

B A b a ir^a-m ° 

matrix The symbol ( | represents the Clebsch-Gordan coeffi- 


cients and f^ (r^) is the form factor depending upon the 
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choice of the neutron bound state wave function and |B n 

1 lsj 

is proportional to the interaction strength D q 2 and the 
spectroscopic factor s^sj 

In the case of (d p) reaction the formula (2 29b) 
reduces to 


do 

ao. 


1 53 


2J B +1 

2J a +1 



Lisj_ 


n IS J 
a DvJ 


(2j+ 1) 


(2 31) 


2 3 

where Hulthen function has been assumed for the deuteron 
wave function Normally only one value of lsj is required 
and therefore (2 31) simplifies to the form 


der 

d^ 


2J A +1 


s „ ls J 
s lsi °DW 

(2 J + 1 ) 


Compound Nucleus Reactions 

A discussion of formation of compound nucleus would 
really need the consideration of the time development of the 
non-stationary state produced m nuclear reactions The state 
at t= - « develops into a bound state at the instant of com- 
pound nucleus formation which subsequently decays providing 
f ] ux m unbound channel m an exponential manner In the 
compound nucleus model it is assumed that the Y k oun a compo- 
nent of the stationary state wave function predominates within 
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the nuclear volume to provide the boundary conditions at the 
nuclear surface and the stationary unbound wave functions 
predominate just outside the nuclear volume The coupling 
between the two is provided by the boundary conditions The 
compound nucleus has a definite spin and parity ir so that 
we have 


+ s a + J A - J c “ % + S b + J B 
Tr a ir A (-l) 10 = Wg = tr c (2 31a) 

It is usual to introduce channel spins S a such that 

S 1 = S a + J A s p = s b + J B (2 31b) 

As a consequence of restrictions (2 31) the values of angular 
momenta m the channels a and 0 are restricted by 



= (J c + s 0 ) 

h 

0 

i 

a 


IB 

= (J c + ) 

l J c - S 8l 

(2 32) 


Further, because of parity conservation the permitted values 
of l a and lg are either all odd or all even The angular 
dependence of the amplatude A a g m channel 0 therefore, has 
the form 

"Snax , , 

A w8 = z a 1 Y 1 

1 ( even or odd ) m m 
m 


(2 33) 
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The cross-section is proportional to the square of the ampli- 
tude and therefore has terms which are of even order The 
dependence on the azimuthal angle vanishes after summation 
over the magnetic quantum numbers and the cross-section can 
be expressed m the form 


do 

d 


L 

max 

X 

L 

even 


P T (cos 


©) 


(2 34) 


after expansion m terms of Legendre polynomials The angular 
distribution/ therefore are symmetric about 90° m the centre 
of mass coordinates If more than one level of the compound 
nucleus of different parities contribute the summation (2 33) 
will be over even as well as odd values of 1 and the angular 
distribution will then no longer be symmetric about 90° CM In 
the case of contribution from two or more levels of the com- 
pound nucleus the angular distribution changes rapidly with the 
bombarding energy since the relative contribution from diffe- 
rent levels changes 

The angular dependence of the amplitude A , 

s m s m 
s s 

from initial channel s m_ to final channel s ' m is given by 

S 

2 a 

(Blatt and Biedenham) 


^s m s ' m 1 
s s 


x (1 s' m^irig, | J c m c ) i 


1 * -1 J, 


- I (21+1)^ (Iso m g |j c m c ) 


1 * 


U 


Is 1 ' s ' 


m. 


(e,«) 


(2 35) 
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where A 


s m s p i is related to amplitude q , 

s s s rn s in f 

s s 


by 


the relation 


}a 




s m s ' m , 
s o 


/0 ,-2i.2 1 q s m s rn , 

(2im ) s s’ 


Eq 2 35 is the detailed form of (2 33) The summation is over 

all quantum numoers except s m s', m The quantity U is 

s s 

the partial reduced amplitude which is independent of angle 
but depends upon the system energy If a single level of the 
compound nucleus is excited, there is no summation over If 

a statistically large number of levels contribute it is 
expected that the interference terms due to their randomness 
may cancel and the cross-section will be relatively insensi- 
tive to 9 Also, angular distribution will be approximately 
symmetric about 90° m centre of mass system The total cross- 
section is obtained by summing over the final channel spins 
and averaging over the initial channel spins m case of un- 
polarized beam 

The calculation of the quantity U is to be carried out 
as a problem of nuclear structure m the region of high exci- 
tation energy This leads to the prediction of the position, 
angular momentum/ parity and other properties of the levels of 
the compound nucleus The compound-nucleus model developed 
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so far does rot attempt to calculate these from first princi- 
ples However, assuming that these levels exist, their effect 
on the nuclear reaction process can be predicted ,71th some 
simplifying assumptions I f at excitation energy E such a 

level exists the probability of finding the incident particle 
a within the nuclear volume is expected to be high since at 
this energy a bound state is predominantly excited The condi- 
tion of high probability within the nuclear volume is satisfied 
when the slope of the radial wave function rR^(r) of the un- 
bound state is zero at the nuclear surface The energy at 
which a level of the compound nucleus is formed is therefore 
also the energy at which the logarithmic derivative f^ at the 
nuclear surface is zero The general quantum mechanical theory 
of scattering m the presence of absorption shows that the 
cross-section for scattering and for compound nucleus forma- 
tion peab for this boundary condition It is therefore possi- 
ble to understand appearance of resonances m terms of levels 
of compound nucleus by introducing boundary conditions at the 
nuclear surface The behaviour of cross-section m the 
neighbourhood of resonance is found to be of typical dispersive 
type as given by the Breit-Wigner formula This feature is a 
consequence of the fact that the logarithmic derivative varies 
reasonably smoothly with energy m the neighbourhood of its 
value zero In that case it can be expanded m power series, 
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about the resonant energy and the retention of the linear term 

alone directly leads to Breit-Wigner formula which is of the 

•p 24 
form 


U ^ 



(E-E n )+ lA 


(2 36) 


The cross-section is proportional to U and for the 
shape of the resonance the parameter A is found to be equal 
to r/2 where r is the full width of the resonance at half 
maximum According to uncertainty principle the width is pro- 
portional to the rate of decay of the compound nucleus The 
rate of decay can be further divided into partial rates of 
decay into various specific channels Ue can thus define 
partial widths F T etc corresponding to these partial rates 

P U 

such that T = S Fg The energy dependence of the quantity U of 
Eq (2 35) is 


J c ^ r nls ^nl's' 

u is I’s' = tt; 

( E " E n } + “ 


(2 37) 


Substituting Eq (2 37) m (2 35), summing over all quantum 
numbers the amplitude q has the form 


q = I 


n 


2n 

{(E-E n ) + 



(2 38) 
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where n represents different resonance energies, q n represents 

dependence on parameters other than energy and the phase factor 

is includes m a 

n 


Contribution of Dire ct and Comp ound Iiucleus Reaction s 

As mentioned earlier the actual problem of finding the 
asymptotic behaviour of the total wave function which includes 
the bound and unbound parts with the coupling between the two 
is treated m terms of models The amplitude is assumed to be 
made of a resonant part and non-resonant part This has basis 
m terms of semi-classical consideration and is applicable 
when reasonably resolved Breat-Wigner type of resonances are 
found riding over general background m the measured cross- 
section This is particularly applicable m the case of reac- 
tions liXe the (d p) reactions on light targets where well 
recognized stripping patterns of angular distributions are 
observed at a bombarding energy of about 4-5 Me/ while at bom- 
barding energy of about 2-MeV the stripping patterns appear m 
a somewhat distorted form and resonances appear m the excita- 
tion curves 

In such cases it is possible to assume that the ampli- 
tude q is given by 


cr(E,9) = 


{(E-E n )+ fa; 


(2 39) 
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TTR 

where the non-resonant amplitude q nay be considered to be 
energy insensitive m magnitude and phase over the v ldth of 
the resonance though it may depend uoon 0 The phase of the 
resonant part of the amplitude varies with energy as deter- 
mined by Breit-Wi Tier formula Due to this energy dependence 
it is possible to empirically separate the non-resonant part 
of the cross-section from the resonant part by analysing the 
experimental excitation curve data at each ancle It is then 
possible to reconstruct the angular distributions of the non- 
resonant mechanism alone and compare these with the available 
predictions of the direct reaction theory 


CHAPTER III 


EXPERIMENT 


Accelerator 

The present work was done with a two-Mev Van-de- Craaf f 
accelerator available at the Indian Institute of Technology 
Kanpur It is a standard two-Mev machine of High Voltage 
Engineering Corporation Burlington Mass USA 

The accelerator was purchased without the gas drying 
system which was designed and assembled locally The system 
provides a fresh charge of nitrogen and carbon dioxide mixture 
with a dew point of -60°F m about ten hours which is then 
transported to the pressure tank at 350 psi to provide insu- 
lation of the high voltage dome A water chilling plant was 
similarly locally assembled 

The ion source is a standard radio frequency source 
which is at present used with hydrogen and deuterium as source 
gases The ions are focussed and accelerated down the beam 
tube to provide total currents upto about 150 & amperes with an 
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energy resolution of about lo KeV The analysing and switch- 
ing magnet for beam handling was under construction while this 
wor 1 was m progress Therefore a laboratory magnet provid- 
ing aoout 10 kilo-gauss of magnetic field was used to separate 
the molecular beam from the atomic beam 

The voltage on the dome was measured by means of a 
generating voltmeter (CVM) fixed inside the pressure tank It 
consists of a rotor which alternately shields and reveals a 
stator to the electric field produced by the charge of the 
dome P n AC signal is therefore generated on the stator which 

is amplified rectified and measured The output DC signal is 
proportional to dome voltage The absolute calibration of the 
beam energy was done by locating the Li (p n) threshold expec- 
ted at 1 88 MeV and the 27 Al(p y) resonance expected at 0 992 
MeV The neutrons were detected with a BF^ counter embedded 
m paraffin and Y-rays were detected with a Nal crystal of 
dimension 5'x4' mounted on photomultiplier The final dome 
voltage was read on a digital voltmeter The average fluctu- 
ation m the dome voltage during the period of data collection 
was +3 KeV 

Target Chamber 

A scattering chamber with a diameter of 10” was designed 
for charged particle experiments and is shown m Figure 1 A 
system of two slits was used to collimate the beam to a narrow 
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spot of two millimeters m diameter and a third slit with a 
larger aperture was placed nearer the target as an anti scatte- 
ring baffle The target was mounted on a brass holder which 
could be raised lowered or rotated without breaking the vacuum 
by adjusting the clamps on the top of the scattering chamber 
The centering of the target could be ensured by visual checks 
of the fluorescence through the perspex window provided for the 
purpose The scattering chamber had the provision of holding 
a number of solid state detectors at any desirable radial dis- 
tance from the target and their angular position could be read 
and varied without breaking the vacuum with a knob at the 
bottom of the chamoer Arrangement for a monitoring detector 
at a fixed angle was also made An electrically lsulated Fara- 
aay cup was provided to collect the transmitted beam This is 
fitted with a secondary electron suppressing ring which was 
kept at -150 volts to prevent leakage of charge through secon- 
dary electron emission from the Faraday cup The output of the 
Faraday cup was fed to an Eloor current integrator which 
measures the current as well as the integrated charge The 
chamber was provided with an evacuating port and a vacuum of 
10 -6 mm was maintained using the evacuating system of the 


machine 
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Detector 

For the detection of charged particles surface barrier 

solid state detectors were used They had a depletion layer 

sufficiently thicl to stop 6 MeV protons and 24 MeV alpha 

particles These were obtained from the Oak Ridge Technical 

Corporation (0R1FC) Oak Ridge Term USA. and some were 

obtained from the Bhabha Atomic Research Centre Trombay, 

Bombay Some detectors made m this laboratory were also used 

for preliminary studies The detector signals were fed to ORTEC 

model 101 and 102 low noise preamplifiers and model 201 and 202 

ORTEC amplifiers The output of the amplifier was analysed 

using a 512-channel multi-channel analyeei (MCA) purchased from 

Nuclear Data Corporation The resolution of the detection sys- 

2 41 

tern was checked with the Am alpha source and was found to 
be 0 QZ for the 5 48 MeV alpha particles It was improved to 
0 7/ when onl^ the central portion of the detectors was used 
with a detector oollimation system A sliqht effect of degra- 
dation of resolution, due to neutron irradiation was observed 
during tne course of (d p) experiment The MCA itself was 
calibrated using a linear pulser provided m the ORTEC ampli- 
fier system 

Target 

12 9 

The targets of C and Be were made by vacuum evapo- 
ration technique The system was assembled m this laboratory 
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and xs shown in Figure 2 The electrode system couplings and 
the liquid nitrogen trap were fabricated locally It was possi- 
ble to obtain a vacuum of lo” ® mm m the evaporation dome and 
a current of about hundred amperes could be used for evapora- 
tion of the material 

For the preparation of carbon targets two carbon elec- 
trodes (99 9/ purity) were fitted to the live posts The tip 
of one was tapered and that of the other was a plane surface 
The tips were brought m touch with a rotary seal arrangement 
to strike an arc ^ tantalum boat was fitted to two other 
live posts and a small quantity of sodium chloride was placed 
m the boat Glass slides etched with hydrofluoric acid and 
cleaned with acetone were held at about 15 cm above the boat 
for deposition A clean glass dome was placed above the whole 
system and the dome was evacuated The sodium chloride was 
first evaporated by passing current through the boat and the 
glass slides were thus coated with thm-layer of sodium chloride 
An"arc was then struck between the carbon electrodes which led 
to the deposition of a carbon film over the sodium chloride 
film The slides were then removed and slowly immersed m dis- 
tilled water which dissolved sodium chloride The carbon films 
floating on the surface of water were then picked up on rectan- 
gular brass mounts with circular holes m them 
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For preparing the beryllium targets carbon films were 
first prepared as packing and the glass slides were replaced 
bv these mounted films for evaporation of beryllium Beryllium 
(99 99 °/ Tunty) v.as evaporated using tungstun filament instead 
of a tantalum boat since beryllium forms an amalgum with tan- 
talum The beryllium targets with carbon backing were used m 
the experiment and the additional peaks m the reaction spec- 
trum arising out of beryllium reactions could be easily detec- 
ted since the carbon peaks were already studied The beryllium 
peaks of interest were well separated from the carbon peaks due 
to the difference j.n the Q values Attempts to prepare a self 
sun-porting beryllium film did not succeed since the beryllium 
flakea-off when attempt was made to peel off the film from the 
glass slides 

The relative measurement of the thickness of the carbon 
targets used was obtained from their respective yields of 
(d p ) reaction The absolute thickness was found by weighing 
a relatively thick film m a microbalance and also by studying 
tne energy loss of 1 90 MeV protons by studying the shift in 
the threshold of Li (p,n) reaction as siown m Figure 3 The 
two methods gave consistent results to within ICP/ The carbon 
targets used had a thickness between 3-16 KeV for 1 9 MeV protons 
In the case of beryllium target the thickness of the backing 
carbon film was first determined ana that of beryllium film 




/ 



Target thickness measurement 


o Without target in the path. 

□ Carbon target in the path. 

* Carbon + Beryllium target in the path 
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was obtained by measuring the energy loss through the two 
films together Beryllium targets used had a thickness of 
3 to 7 KeV for 1 9 MeV protons To improve an overall accu— 
racy, thick-target thickness was determined 


Measu r em ent s 

For alignment of the beam, the Faraday cup was replaced by 
a pyrex disc and the position of the scattering chamber was 
adjusted to ensure that the beam spot was at the centre of the 
pyrex disc The target was then lowered and the passage of the 
beam through its centre was ensured by checking visually the 
bright spot on the target appearing due to beam bombardment 
The pyrex oisc was then replaced by the Faraday cup The zero 
degree position for the detector was checked by observing the 
yield of the oro tons scattered from carbon 12 at the same posi- 
tive and negative angles for several angular positions 


A typical spectrum at 1600 KeV deuteron energy and 
9 lab = - 1 -- 1 - 30 for beryllium on carbon oacKing is shown m Figure 4 
All resolved peaks m the spectrum could be identified from the 
known Q values for the respective reactions For reaction angles 
smaller than 70° lab the elastic scattering cross-section in- 


creased rapidly as compared to reaction cross-sections This 
increased the dead time of the equipment to an undesirable 
extent Therefore a nickel film was placed m front of the 







Channel No 
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detector to just stop the elastically scattered deuterons but 
not the energetic protons from (d p) reactions With nickel 
foil m front of the detector data upto 15°Xjab could be taken 
m the forward direction \ typical spectrum with nickel foil 
m front of the detector is shown m Figure 5 To identify 
peaks m this case loss of energy m Ni foil for each reaction 
product had to be calculated separately 

12 

The excitation function for c(d,p o ) reaction was 

obtained at nine laboratory angles m steps of 10 Ke\7 m the 

dcutcron energy range of 800 to 1850 keV shown m Figures 

12 

6-9 The excitation c-rve for C(d/d) channel, m the same 
range is shown in Figure 10 Further for consistency check, 
angular distributions m the angular range 15° to 160° labo- 
latory angles m steps of 10° were taken at seventeen values 
of bombarding energy in the range 800 to 1800 KeV These are 
shown m Figures 13 to 13 In the low energy range of 500-800 KeV 
no resonances were observed Excitation curve m this range at 
the three laboratory argles and angular distributions at two 
bombarding energies are shown m Figures lo ano 

On ^8c target aata on (d p ) (d t) (d,d Q ) and (d a^) 
could be ob rained Since these reactions have a higher Q value 

than the ^ 2 c(d o ) reaction, the corresponding peaks could be 
studied mspite of carbon backing at most of the angles In 
the energy range of 800 to 1800 KeV excitation functions for 
9 Be{d p ) reactions were taken at nine angles These are 
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shown in Figures 6 , 15 and 16 For consistency check angular 
distribution m the angular range 15° to 160° laboratory 
angles m steps of 10° were taken at eighteen values of bom- 
barding energy m the range 500 to 2000 KeV as shown m 
Figures 17 to 20 In addition excitation furctions for 9 Be 
(d p l > reaction m the range 800 to 1800 KeV and excitation 
function for ^Be (d ) m the same bombarding energy range 
are shown m Figures 21 and 22 Excitation curve for 9 Be (d t ) 
m the range 800—1800 KeV is shovn m Figure 23 

It is estimated that the cross-sections have been mea- 
sured with an energy resolution of 0 5/ The relative error 
m the cross-sections is about 2/ and the error m absolute 
cross-section is about 12/ m carbon and 15/ for beryllium 
In the present analysis the accuracy m relative cross-section 
is of primary importance 
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CHAPTER IV 


ANALYSIS OF DATA 

The excitation functions of 12 C (d,p ) reaction leadi n 

13 ^ 

to the ground state of c at nine angles is shown in Figures 

6 to 9. The excitation function of l2 C (d,d) reaction at two 

backward angles is shown in Figure 10. The striking feature 

in this data is the appearance of correlated peaks at all 

angles in (d,p ) reaction as well as in the (d,d) channel, j. 

example there is a rather well separated peak appearing at 

about 1450 KeV in all these excitation curves. Interference 

effects are obviously present and these lead to the difference 

of details of the excitation curves at different angles. Thu«=> 

the 1450 KeV resonance is sharpest at 140° in the (d/d) c h?,r„ 

and is weak at 80° in the (d # p) channel. Other peaks which 

less marked stand out to different extent at different ang]_ es 

Thus the peaks in the vicinity of 930 KeV, 1310 KeV, 177c p 

and a broad hump around 12 00 KeV can be noticed, which is 

sis tent with the previous data. The data indicates the 
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presence of a resonant orocess with a level density which is 
not too hiwh ^ms is expected for the relatively light com- 
pound nucleus N at an e citation eneigv of about 11 5 MeV 
The data obtained m the past, has been analysed by assuming 
the existence of compound nucleus levels at appropriate exci- 
tation energies 

The anqular distribution of the ^^c(d,p Q ) reactaon m 
the energy Lange covered is shown in Figures 11 to 13 The 
noticeable feature of most of these angular distnburions is a 
maximum at about 35° and a strong backward peaking The maxi- 
mum at about 35° is typical of the direct deuteron stripping 
reaction while the backward peaking is not The angular dis- 
tributions on and off 1450 FeV resonance are shown m Figure ?4 

The angular distributions do not change much m their general 

' | 

shape, however, the backward cross-section is large at the pea^ 
energy The dashed curves shown are tnose predicted on the 
basis of P] ane Wave Butler theory with l n = 1 for the captured 
neutron The fits obtained for the 35° maximum are satisfac- 
tory However the theory grossly underestimates the relative 
cross-section at backward angles and at the minimum of tne 
angular distribution The available data at higher energies 
shows that the backwaid peaking m the angular distribution 
decreases until more familiar stripping patterns are observed 16, 18 

(shown m Figure 25) 
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The excitation function of "'Be(d p ) reaction leading 
to the ground state of ^°Be at nine -angles is shown m Figures 
6/ 15 and 16 The excitation function of nrotons from the 
'Be(d p^) reaction leading to the first excited state of ^°Be 
is shown m Figure 21 The excitation functions of Q Be(d ct ) 
and Be(a t) leading to the ground state of Li and b Be respec- 
tively are shown m Figures 2 2 and 23 Correlated weah reso- 
nances or antiresonances are observed m the excitation curves 
mentioned above at deuteron energies of 950 keV 1080 KeV 
1200 KeV 1300 KeV, 1450 'eV 1620 KeV and 1760 KeV The reso- 
nances m all these cases are weak and show that compound nucleus 
effect is small as compared to those in the case of C (d,p ) 
reaction The excitation energy of the compound nucleus 

q i 14- 

-'Be + d = B is about 17 KeV which is higher than in the M 

compound state Many channels are open m this case for the 

decay of the compound nucleus and therefore contribution to 

any one channel as expected to be small On tne other hand 

stripping process if present will contribute predominantly 

19a 

to (d, p) and (d n) channels Earlier Canavan has analysed 
^Be (d, p^) data m this energy range assuming that compound 
nucleus process alone is responsible for the reaction mecha- 
nism Later groups concluded that the reaction was taking 

place mainly through direct mechanism and tried to fit data 

„ 20 

using the stripping theory 
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The angular distributions of the Be(d p ) reaction 
n t about 2 HeV is shown m Figures 17 to 20 In the forward 
angles a maximum at aoout 35° is observed at most of the ener- 
gies Also an equally strong bad ward peaking throughout is 
observed The pea" 1 at 35° is typical of stripping reaction 
mechanism The angular distribution changes a little as the 
energy is changed Fowever, the relative forward to bad tfard 
cross-section increases with the erergy A typical angular 
distribution at 1310 FeV 1445 FeV and 1620 KeV is shown m 
Figure 2 6 Dashed curves are Plane Wave Butler curves with the 
neutron captured m the ' p ' orbit The fits are good in the 
forward angle but predictions of relative cross-section remains 

poor m the backward angles The angular distribution data at 
2 S 2 6 

higher energies (Figure 27) shows considerable deciease 

m the relative backward cross-section ana a pure stripping 
t^pe of pattern with a small backward cross-section is observed 
at these energies 

In view of these observations it was felt that the 
possibility of coherent contribution from both the resonant 
and non-resonant processes should be investigated although as 
was done m the past, it is possible to fit the angular distri- 
butions with Legendre polynomials and explain them in terms of 
contribution from a number of compound nucleus levels exclu- 
ding apriori the possible non— resonant contribution 
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In this connection the ancrular distributions obtained 
\\pre fitted oj us m terns of Legendre polynomials It was 
found thrnt a satisfactory value of ^ m the least square 
fitting could be obtained by using both even and odd terms 

unto a maximum order of l given by L = 6 m the case of 

X P Q 

C (d,p ) and L = 7 m the case of Be(d p ) We have 
o max o 

L <21 
max — 

x/here 1 is the smaller of the two values l a and 1^ In the 

12 9 

case of C we therefofe have l a _> 3 and m the case of Be 

l a >_ 4 To explain the observed distributions on the basis 
of the compound nucleus theory rather 1 iga partial waves m 
tne incident beam have to be mvohed On the otner hand it 
is e oected that such high partial waves would lead to stri- 
pping mode rather than compound nucleus formation 


J 


Jt 


Sepa rati on of v" 1 J esonan t_ Amplitude 

The separation of compound nucleus part from excitation 

1 6 

curve has been earner suggested by Bonner et al and Yoshida 

The separation of non-resonant part from the excitation curves 

is done on the basis of its smooth energy dependence The 

differential cross-section m terms of the amplitudes q sm g m 

s s ' 

i<- given by the relation 


do 

an 


s m s 1 in f 
s s 



I q , | 

1 bin s m , 


2 


(4 1) 





79 


The amplitudes q , is split into a non-reso- 

s n s rr * 
s s 

nano (NR) and a resonant (R) part \t a gi'/en angle 9 we 
have 

q(E a ) = Q TP (E a ) 4 - q R (E a ) (4 2) 


where the channel quantum numbers have been dropped for the 
convenience of writing the expressions The resonant amplx 
tude has the following energy dependence 


q ? (E a ) = l 


r> 3-6 

C R e n 
n 


n (E a" E n )+ 1 V 2 


(4 3) 


where the summation is over all tne contributing resonances 

P 1 6 

and C R is independent ox energy E a The quantity e is the 

th 

relative phase due to n resonance 


Uhen onli one level of tne compound nucleus is consi- 
dcr 1 and m Eq (43) contribution from single resonance is 
considered the cross-section is given by 


dcr 

df 2 


NR 

q + 


° C R 
o 


s m s m , 
s s 


K, 


(E -E )+ ir /2 
a o o 


2 


(4 4) 


The amplitude q 1P is assumed to be constant over the width of 
the resonance. Equation (4 4) therefore becomes 


t 
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do 

dP 


sms m 
s s 



! cr 


HR 


2 



‘W + 



/ 4 


+ | q FR C R | 2 COS j o (E a" E o )+ 2 F n4 

(F a -E o ) ? + r Q 2 / 4 

. x k 2 ( 4 5 > 

U ° re factor ~ has been absorbed m the amplitudes To 
evaluate the cross-section ~ the Eq (4 5) lg sunned over 

the f nal channels and averaged over the mitxal channels 
The relative contribution of each term in Ecr (4 5 ) does not 
change by this averaging Carrying out proper summation and 
re-defining the constants 3 n (4 5) after terms of similar 
dependence are grouped together we get 


do 

dP 





<^-s o ) 2 + r o 2 /4 


+ 


3* ( (%-c o ) 2 +r o 2 /4) 


(4 6 ) 


The constants D C Q and are independent of energy 
The first term gives the value of non— resonant part 

ox the cross-section The Eg (4 6 ) can be re-written as 


Y = F P , + P 2 f 2 (E E o r o> + P 3 f 3 (E E o r o’ 

(4 7) 

The Eq (47) becomes a basis for separation of non-resonant 
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contribution from total cross-section, when contribution from 
cnl/ one resoranc" a.s considered The constcints and 

are treated as parameters lxi least squares fitting of this 
expression 7) 


If cortr bution from two levels of the compound nucleus 
is considered Err (4 3) is re-wriut Q n as 


q(E a ) 


q t,p <V + C ° e 


R _. 15 0 


P 16 1 
C 1 e 


(E a -E o )+i r o /S (E^ep + ir/2 


(4 8] 


The cross-section ,-x . is given by 

a« 's m s m 
s s 


d£_ | __ 1 

dP 's m s 1 m T , 2 
s s K a 


l 6., l (S-, 

^ r ^ 1 l 

Q NR (E a) + 1 + 1 


(E a -E Q )+ i%/2 (E a -E 1 )+ ~\/2 


*” NR i 2 

q I + 


K 


R 


lh R f 


a 


(E a" E o )2+r o 2/4 


. HR R 

+ c] Co X 


{ 2 cos 6 ^ ( E a -E i' ) 

2cos S Q (Ij, -E a ) + 2 p/2 Bin 6^ HR „ R + 2 V 2 sln ! 1 5 
( — ')+ a r-i 


(e„-e^) 2 + r 2 /4 
' a o o 


1 (E a -E 1 ) 2 +r 1 2 /A 


+ 2C 


R 

o 



(COS !s 0 -«,_)((E a -E o ) (E«-E 1 ) + r 1 r o /4) 

-sm ( « n - «, ) ( p/2 (E„ -E 0 ) -r o / 2 < E C.- E 1 > ) } ; 

[■ — J 

(<e 0 - e o ,2+ o 2/4><(E «' E 1 )2+ V 4> 


(4 9) 
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Eq (4 9) is further simplified by collecting terms of ssular 
gy dependence After summing over the final channels and 


averaging over the initial channels the differential cross 
section is given by 


da _ D 
dft ~ E_ 


+ 


E et ( ( E a“ E 0 ) 2 + r Q 2 /4) (E a -E ) 2 +r 2 / 4 


+ 


o' "o' E a ( ^ a “ E l^ r i 2 A) 


+ 


(E a El ) 2 + T^/4 E a ( (E a _E o ) ? + T q 2 /4)( (E a — ) 2 + 


C 5 

^ E a“ E 0 )2 + r Q 2 / 4 ) C(E a ~ e^) 2 + r^/ 4) 


C (E a -E 0 )2 + r o 2 /4) C(E a -E 1 ) 2 + r ± 2 /4) 


(4 lo) 


Ihe corstn ts' n crr'cr'r-'^.-i_ , 

'o l ^2 j C 4 E 5 e 1 ’ iave been re-defined 

t-o include all other constants The relative phase 6 ] _ and & 2 
are also treated as constants for a given © and meraed m these 
constants The function (4 10) can be rewritten as 


D 


Y E* + C o f o (E a E 0 V + C l f l^ E a V r Q ) + T ± ) 

+ C 3 f 3^ E a E 1 V l ) + C 4 f 4 ( '\ E 0 /E l r 0 ^ 

+ C 5 f 5 (E a ,E Q E i r Q T ± ) + C 6 f 6 (E a E q E ;[ T q 


(4 11 ) 
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ixpz ess ions (4 7) and 11) could now be used for least 
squares fitting to the experimental excitation curves taben 
at- various angles 

While in nrmciple it was possible to fit the excita- 
tion curves m the entire energy reaion stucied by considering 
t^e contrioution of each resonance everywhere, the number of 
terms m the fitting function and the number of parameters 
become large Instead , a few regions of the excitation curve 
around relatively isolated resonances were studied In the 
first step the contributions from the resonance under consi- 
deration and the contribution from the non-resonant process 
were considered The formula (4 7) was therefore used for 
least sauares fitting m this case T ext the contribution from 
one nciqhoourmg resonance was added and the least square 
fitting was done thi £ tine with formula (4 H) This was 
extended to consider the contribution to a single resonance 

by th^ non iLsondi process and by the two neighbouring reso- 

do 

nance-. The a~m at each step was to obtain the value of — 

du 

given by the first, term of the fitting functions at each 
angle in the cases studied it was found that the non-reso- 
nant cross-section as obtained by using function (4 7) changed 
when function (4 11) was used instead but changed little 
thereafter Therefore for theextraction of — it was 
found sufficient to consider the contribution to the mam 
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resonance from the non-resonanf process and one neighbouring 
resonance The resonant energies and the rfi^ths of resonances 
wore varied along ith the coefficients j.r ■'he fitting func- 
tion to obtain the best fit By this nethoa it was found 
that the excitation function m the neighbourhood of 1450 Fev" 
m the case of C (d p ) could be accounted for by resonance 
of width 10 I eV located at 1445 KeV, a non-resonant contribu- 
tion and the contribution from a resonance of width 40 keV 
located at 1310 The fits obtained m this region of exci- 

tation curve are shown m Figure 28 The non-resonant contribu- 
tion at 1445 Re/ v,as determined at each angle at wmch excita- 
tion curves were obtained It was then possible to plot the 
angular distribution of the non-resonant part alone at 1445 KeV 
bombaidmg er rgy This extracted non-resonant cross-section 
was then cornua 1 " ed with the DWBA theory of stripping This 
nroceuULe wa r repeated at the 930 Ke\ resonance whose width 
was determined as 110 KeV a contribution from neighbour reso- 
nance located at 1160 KeV was considered Wxdth of this reso- 
nance /as deter med to be 190 KeV The non-resonant cross- 
section was a 1 so extracted at 1770 KeV with 140 KeV width 
The extracted angular distributions due to non-resonant process 
are shown m Figure 29 The non— resonant contribution accounts 
gbout 20 °/ of the observed cross— section 
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9 

In Be the non-resonant cross-section contribution was 
determined at 920 KeV 1310 FeV 1450 FeV 1620 keV and 1770 
FeV where widths of resonances were identified as IOO KeV 
50 KeV 80 KeV 50 KeV and 120 FeV The angular distributions 
due to non-resonant contributions are shown m Figure 30 The 
non-resonant cross-section m this case accounts for about 
9 oy of the observed cross-section 

DWB4 Calculations 

Ihe extracted angular distribution due to the non-reso- 
nant process was compared witn the direct reaction theory of 
the stringing orocess The DWBA calculations were carried 
out with program DWIJCK ana the computer BESh-16 at Bhabha 
Atomic Research Centre Trombay Bombay 

The dis orted wave program DWUCK calculates the diffe- 
rential cites s-section as a function of angle 9 for a given 

Is 1 

energy of the incident particle The quantity cr DW as 
given bv Eg (2 29b) is calculated for a given value of 
orbital angular momentum transfer 1 sum transfer s and 
total angular momentum transfer j The distroted wave X L 
tor partial wave L is calculated by solving the following 
differential equation numerically 

+ k ? - L - ^ ~ + ^ (V B )] X L (r> =0 (4 12) 

dr r Tl 


y 

h 
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f 






w 1 ere 

v f = u c + u so 

and the potential U^(r)i S aiven by 

d f(X_) 

U N (r) = 7 P f ' 7 p) + 1W X (4 13) 


where W R is 
nary part 


real part of the potential and W is the imagi- 
The function f(x) is given by 

f (X) = (1 + e 7 ) 
r-R A 1 / 3 


x = (—2 ) 

a = Diffuseness parameter 
A = Atomic Number of the nucleus 


The potential U gives Coulomb interaction due to charged 
sphere of radius R and charge Z^ and is given by 


2 A e 

U = - - - -- r > R 

c r c 


Z, Z A e r 2 

(3 ' ^ r — 
C 


(4 14) 


The interaction U c _ is proportional to spin orbit interaction 

oU 


L S and is given by 


U S0 " W R 45 2 l 


VSOR 1 d f(X R } 


dr 


L S 


(4 15) 
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where VSOR comes from the Thomas Spin Orbit term and has a 
value of 25 for j rotons and neutroi s The wa^e function 
X T (r) satisfies th^ boundary condition m the asymptotic 
region where — ? 0 and also inside the nucleus where 

X. (O) — > 0 The ^(r) are computed by numerically integra- 
ting the differential ecruition 12) and thus values of 

X (r) for a set of discrete values of r is obtained m the 
entrance and the exit channel with proper choice of potential 
parameters m the channel concerned In the case of stripping 
reaction where a particle is captured m the outer orbit, the 
form factor is obtained by calculating the bound state wave 
function for the captured particle such that it reproduces 
correct binding energy Thus for stripping 


= K Nl (r) 

where R^Cr) is the normalized radial wave-function of the 
captured particle After the calculation of form factor 
f (r), and distorted wave x T integral xn Eq (2 30) is 
obtained by doing numerical integration Finally program 
calculates the quantity Jo# 1SJ as given bv Eq (2 29b) by 
inserting proper values of Clebsch Cordon coefficients and 

spins of the system 

There is more than one choice for the form of potential 
U(r) m Eq (4 12) 13a . Che commonly used form is given in 


r 
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Eq (4 13) Here real potential is Wood-Saxor type and imagi- 
nary potential is derivative of this function 

The calculation of integral m Fq (2 30) is done using 
zero-range approximation The DWUCK program also has facility 
for taking into account correction due to finite range effect 
and is calculated m the following way 

The overlap of wave function ^ and <j> a is given by 


MX) = 1 d t o » b * * a 


and the founer transform G(K ), is defined as 


G(K 1 2 ) 


lk X 

/ e D (X) dX 


The zero range approximation is given by 


D q = G(0) 


and finite range correction parameter FNRNG is evaluated by 

2 2 

taking first term m the expansion of G(K ) about K — O, that 


R 


2 


1 

G(K 2 ) 


9G(K 2 ) 

MK 2 ) 



Thus m the reaction where 

a = b + c and B = A + c 
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the finite range correction 


-ls computed by the quantity 


W(r) = { 1 + JL. 


n 




M. 


R 2 (Eh" v (r) + E - V (r) - E -V. (r))> 

^ o a a 


(4 16) 


Potential U N (r) m differential Eg (4 12) is a local 
potential The program DWUCK ta] es into account correction 
due to non-locali cy The effect of the non-locality is taken 
into account by multiplying a wave function for an equivalent 
local potential U N (r)by the factor 28 


f (r) 


(1 


1 

A 2M 

4 ^2 N 


(r)) h 


(4 17) 


ln the channel under consideration Here S is range of 
ron— locality and is 0 85 for protons, 0 54 for deuteron and 
0 2 to 0 3 for 3 He 

The effect of finite-range correction and non-locality 
m general results m reduction of the cross-section but the 
over all change m the shape of angular distribution is not 
expected 

In calculating the DWBA cross-sections wi th DWUCK, the 
principal input consists of the optical model parameters in 
both the incoming and outgoing channels The values of optical 


parameters which provides satisfactory fit at high energy of 
25 29 3 0 

about 5 MeV were used as mtial values in fitting the 

ectracted curve These values are given in Table 1 


Target 


Incident 

Particle 


Table 1 


Real Potential 

__ — 


Imaginary Potential 
( Derivative type) 


1 25 0 47 


1 

41 

O 

65 

1 

25 

0 

47 


In the first place, since the deuteron bombarding 
energy is lower than the Coulomb barrier it was thought that 
the nuclear interaction of the optical potential may not be 
necessary and Coulomb distortions alone would be sufficient 
for the calculations In this case it may he noted that 
according to Eq (2 18) the fully distorted waves including 
nuclear di stortions,must be used m the final channel 6 
although the wave function 'F + could be tried as the incident 
wave with the Coulomb distortion alone This is necessary if 
V 0 -Vo is to be equated to V Preliminary DWBA calculations 
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with only the Coulomb interaction m the incident channel 
were carried out The calculated angular distributions were 
somewhat similar to the Butler curves with higher cross-sec- 
tions predicted at backward angles One such curve for 
S 

Be (d p Q ) at 1445 KeV is shown m Figure 31 The fit for 
the backward angles is poor * more serious discrepancy is 
however m the absolute cross-section "here the predicted 
value is about 15 times the observed value It was therefore 
concluded that a nuclear distortion of the incident wave can 
not be neglected 

In the DWB\ calculations for 12 C (d p ) the optical 
parameters were varied about their initial values Further 
calculations were made with and w thout spin orbit term with 
finite range correctior and possible correction due to non- 
local effeccs were also investigated The inclusion of spin- 
orb t term m che optical potential increased the predicted 
bad angle croso-section The inclusion of corrections due 
to finite rmqe and non-local effect did not alter the shape 
of the angular distribution significantly but reduced the 

Q 

overall cross-section (Figure 32) For Be(d,p Q ) reactions 

only finite range correction was included Figures 33 to 36 

show D rT BA nts to C (d p ) with onlv spin orbit term Figures 

37 to 41 show DWBA. nts to ^Be(d p ) The optical Parameters 

12 9 

which gave the best fit m the case of C and Be are given 
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m Table 2 


Table 2 


Target 

Incident 

Particle 

Real 

Potential 

Imaginary Potential 
(Derivative type) 

W R 

r R 

a R 

W I 

r l 

o,_, 

I 

Q 

"Be 


D 

-160 

0 9 

0 9 

-12 

2 1 

0 ^7 

r = 
c 

1 3 

P 

-49 

5 1 25 

0 65 

-7 0 

1 25 

0 65 

1? c 


D 

-69 

1 6 

0 61 

-4 0 

1 41 

0 65 

r c = 

1 3 

P 

-52 

4 1 45 

0 65 

-5 4 

1 22 

0 65 


It may be noted that the optical Parameters given m 
Table 2 are close to those given m Table 1 for higher energy 
data 

Expression (2 31) for the direct cross-section for 

Is 1 

(d p ) reaction involves the spectroscopic factor s 
A comparison of the DW3A calculation with the experimental 
results provides a value of spectroscopic factor s which 
xs to be compared with that obtained frcm nuclear structure 
calculations ^ The values obtained frcm the P/fBA analysis 
of the data and those obtained from nuclear structure calcula - 
tions are shown m Table 3 
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Table 3 





Experimental 


Caicd 

Energy (KeV) 

920 

1310 

1445 1620 

1770 

9 Be(d p o ) 

2 2 

1 9 

2 0 19 

1 7 

2 36 

12 C (d p Q ) 

3 0 

1 2 

1 1 

1 2 

0 613 


It may be noted that the DWBA analysis of the unextrac- 
ted experimental curve would have led to a spectroscopic factor 
which would be larger than the theoretically expected values 
by a factor of about 8 for the case of C(d,p Q ) 

We may summarise the results of the analysis as follows 
The excitatior functions m 12 C (d p Q ) reaction show resonances 
which appear at various angles and also m the (d d) channel 
The angular distributions indicate the possible existence of 
a non-resonant process also It is possible to extract the 
non-resonant cross-section from the observed excitation curve 
on the basis of known energy variations of the resonant ampli- 
tude (Breit-Wigner formula) and the slow energy variation of 
the non-resonant amplitude Ihe non-resonant contribution is 
about 2 0 0 / The extracted angular distributions have been 


ompared with the DWBA. calculations made with optical para- 
eters close to those used at higher energies The agreement 
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xn the shape and absolute value is satisfactory The results 
Be(d P Q ) are similar Tne resonances are however weaber 
and the direct contribution is found to be about 90? This 
is understanc able since due to higher excitation of the com- 
pound nucleus 1] b , more chcmnels are open and a corresponding 
smaller percentage of the compound nucleus contribution is in 
the proton channel Although the fit m the backward direction 
m the case of ''Be is not satisfactory the forward cross- 
section m tne angular dependence and magnitude is well 
reproduced 
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CONCLUSION 


The air, of this worh was to understand the reaction 
mechanism of <d,p> reactions on light targets at sub-Coulomb 
bombarding energies in the range of 0 8-2 0 MeV The targets 
used were 9 Be and l2 C The features of the stnpoing mecha- 
msm such as the small angle peak for 1=1 capture which are 
easily identified at about 4 MeV bombarding energies were 
still present thou ah not as distinct as m the high energy 
data indicating the possible presence of the stripping mode 
at low energies The excitation functions however showed 
resonances typical of compound-nucleus-formation mode Assuming 
a coherent resonant and non- resonant contribution it was possi- 
ble to separate the non-resonant cross-section and compare it 
wibh that predicted on the basis of DWBA calculations of the 
stripping mode Prom the analysis the following conclusions 


could be drawn 
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(1) The aor roach of accounting for the asymptotic 
behaviour of the unbound wave-functions , which are coupled 
to the bound wave functions in terms of coherent contribu- 
tions obtained from direct reaction model and compound -nucleus 
model (based on physical considerations) leads to satisfactory 
agreement with the data in the reactions studied 

(2) The contribution of direct reactions to the cross- 
1 2 

section m C (dp) reaction is about 2CP/ and in 9 Be (d p ) 

o 1 o 

it is about 9 CP/ m the erergy range studied 

(3) Even at sub-Coulomb bombarding energies used it 
was founa that the distortion of the incident wave is impor- 
tant and the Coulomb distortion alone does not lead to satis- 
factory results 

(4) At higher bombarding energies of about 4 MeV the 
stripping mode predominates and the optical model parameters 
for the DWBA calculations have been well determined Essen- 
tially the same parameters can be used to fit the angular dis- 
tributions at the energies studied provided that the resonant 
part is extracted out 

The effect of finite-range correction and non-local 
12 

correction m the C (d,p ) reaction is to reduce the cross- 
section m magnitude by about 1 C?/ without much change m the 
shape of the angular distribution 
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( ) Resonances were located in the compound nucleus 
1 4 

N at excitation energies consistent with earlier results 17a 
The widths of the resonances which we obtained are however 
somewhat on the lower side m most cases This difference 
is to oe expected since we have included non-iesonant contri- 
bution m the analysis The spin ana parity assignment made 
previously by neglecting the non-resonant contribution needs 
re- examination 
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